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Abstract 

Pure three-qubit states have five algebraically independent and one algebraically dependent 
polynomial invariants under local unitary transformations and an arbitrary entanglement measure 
is a function of these six invariants. It is shown that if the reduced density operator of a some 
qubit is a multiple of the unit operator, than the geometric entanglement measure of the pure 
three-qubit state is absolutely independent of the polynomial invariants and is a constant for such 
tripartite states. Hence a one-particle completely mixed state is a critical point for the geometric 
measure of entanglement. 
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I. INTRODUCTION 



The concept of quantum entanglement has played an important role in the development 
of modern quantum physics. The entanglement theory has its roots in the key discoveries: 
quantum cryptography [1, 2], dense coding [3, 4], perfect teleportation [5, 6] and quantum 
factoring algorithm [7, 8]. These effects are based on entanglement and all of them have 
been demonstrated in pioneering experiments. 

However, the phrase "based on entanglement" does not reveal the essence of those effects. 
For instance, one would like to know what are the states that can be used as a quantum 
channel for perfect teleportation and dense coding? And why those states are capable, 
while others are not? In the case of bipartite systems these questions have a concise answer, 
namely, the maximally entangled states can perform the task and others cannot [9]. But 
the situation is changed drastically in a multipartite setting. For instance, in the case of 
three-qubit systems among the states applicable for teleportation are: the Greenberger- 
Horne-Zeilinger(GHZ) state [10] which is maximally entangled and a W-class state [11] 
which is not maximally entangled since the only maximally entangled state within W-class 
states [12] is the W state [13]. Therefore in the case of multipartite systems the property of 
being maximally entangled is unrelated to the ability to perform a certain task. 

In the meantime the aforementioned three-qubit states, as well as two-qubit maximally 
entangled states, that are quantum channels for perfect teleportation have two common 
properties. First, in each state the reduced density operator of a some qubit is a scalar 
multiple of the unit operator. Second, all those state have the same maximal product 
overlap [14-16]. Recall that the maximal product overlap g of a pure state \ip) is given by 



where the maximization runs over all product states and the normalization of the local states 
\qM){M = A, B, C) is understood. 

Then one makes a guess that these two properties are interrelated and the reasons of 
the interrelation should be analyzed. Moreover, if a pure state has a bipartite entanglement 
that does not depend on the bipartition and is maximal for all possible bipartitions then it 
should have both these properties [17]. 

In this article we consider pure three-qubit states and prove that if the density matrix 
of one qubit is a completely mixed state, then the maximal product overlap g of the state 
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is l/y/2. The basic point of the proof is the existence of the specific point in the space 
of entanglement parameters [18]. In the case of pure three-qubit states there are five al- 
gebraically independent and one algebraically dependent polynomial invariants under local 
unitary transformations [19]. They can be regarded as coordinates on the space of entan- 
glement types and thus any entanglement measure should be a function of these invariants. 
The three polynomials quartic in state function are the squares of the lengths of the Bloch 
vectors and they play a crucial role. If one of these three quartic invariants vanishes, then 
the sextic polynomial vanishes too. Given that all invariants are independent, the vanishing 
of the sextic polynomial indicates the existence of a specific point. Hereafter these points 
are referred to as critical points. 

The term critical point has a mathematical and a physical justifications. The mathe- 
matical justification is the following. The geometric entanglement measure of quadrilateral 
three-qubit [15] and general W [20] states have been computed analytically and the answers 
show the the gradient of the measure has a jump at these points. Hence they are critical 
points. The physical justification is that at the edge of the region of possible values of a 
quartic polynomial the entanglement of the state is absolutely independent of the remaining 
five invariants and then the state acquires an ability to be a quantum channel for the perfect 
teleportation and dense coding. 

This article is organized as follows. In Sec. II, we review local unitary invariants of pure 
three-qubit states. In Sec. Ill, we formulate and prove the main theorem. In Sec. IV, we 
show the for four- and higher-qubit states the theorem ceases to be true. In Sec. V, we 
discuss our results. 

II. LOCAL UNITARY INVARIANTS OF PURE THREE-QUBIT STATES 

In this section we briefly review local unitary invariants of pure three-qubit states [19] 
and slightly modify original definitions for convenience. Consider a pure state ip of the three 
qubits A, B and C. 

There is one independent invariant quadratic in if; which is just the norm of the three- 
party state and therefore has no physical significance. We set it equal to 1. 

There are three independent quartic invariants, namely Tr(p^), Tr(p|) and Tr(p^), where 
PAiPb and pc are the one-particle density operators of the qubits A, B and C, respectively. 
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We redefine these invariants as follows 

b\ = 2Tr(pi) - (Trp A ) 2 , b 2 = 2Tr(p|) - (Trp B ) 2 , b 2 = 2Tr(p 2 ) - (Trp c ) 2 . (2) 

Of course, the substraction of the constant (T^a) 2 = (Trpe) 2 = (Trpc) 2 = 1 does not change 
anything and the new invariants are independent too. The advantage of the redefinition is 
that the new invariants b A , &b and be are the lengths of the Bloch vectors b A , b# and bp 
of the qubits A, B and C, respectively. Therefore 6a/ = (M = A, B, C) means that pu is 
a completely mixed state and bu = 1 means that pu is a pure state. 
There is one sextic invariant given by 

t = 3Tr[p AB (p A ® Pb)} ~ Tr(pi) - Tr(p B ) - ^(Trp c ) 3 . (3) 

Again our definition slightly differs from the original one since we subtracted the last 
term(equal to 1/4) that does not exist in the original definition. The meaning of the sub- 
straction of the last term is the following. Let us define a bipartite correlation matrix G by 
formula 

dj = Tr(p AB ffi®^), (4) 

where a^s are Pauli matrices. Then t is expressed via the Bloch vectors and the correlation 
matrix as follows 

t=^b A -(Gb B ), (5) 

which does not contain an additional term 1/4. 

The invariant of degree 8 is (up to a numerical factor) the square of the three-tangle 
r [21]. We use it as is. 

The last invariant % is discrete in the sense that it can have at most two different values 
when the other invariants are fixed [23]. It appears as follows. The polynomial invariants 
constructed in Ref. [19] do not distinguish a state and its complex conjugate and therefore 
one additional invariant is needed to specify uniquely a pure three qubit state. However, all 
amplitudes of a three-qubit pure state can be chosen positive when a Bloch vector vanishes 
as will be shown below and then the additional invariant does not play any role here. This 
point will be explained more clearly below. 

We listed all independent invariants that are coordinates on the space of orbits of the 
group of local transformations. An arbitrary entanglement measure, and among them the 
geometric measure, is a function of these variables [19, 22]. 
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III. MAIN THEOREM 



It is curious that there are specific points in the space of entanglement parameters 
{pAi bs, be, t, t, i). The following theorem clarifies the specific points and the physical mean- 
ing of those points. 

Theorem. 

If b A b B b c = 0, then g 2 = ^. 

The main objective of the article is this theorem. It states that if any of invariants b A , bs 
or be vanishes, then the geometric measure given by 

E g {tff) = -2hxg(^) (6) 

is absolutely independent of the remaining entanglement parameters and is a constant for 
these states. 

To prove the theorem we use the generalized Schmidt decomposition that is closely related 
to the geometric measure [24] and hence it is more appropriate to our analyze here. It states 
that a pure state has the canonical form [13, 24] 

\*l>) = a|011) + 6|101) + c|110) + d|000) + e i7 /i|lll), (7) 

where labels within kets refer to the qubits A, B and C in that order. All the coefficients 
a, b, c, d, h in (7) are positive and the gauge phase 7 ranges from — n/2 to ir/2. 

We need to compute explicitly the Bloch vectors and correlation matrix for the following 
analysis. The computation is straightforward and yields 

= (2/iacos7, 2hasmj, d 2 + a 2 — b 2 — c 2 — /i 2 ) , (8) 
b# = (2Kb cos 7, 2/i6sin7, d 2 + b 2 — a 2 — c 2 — /i 2 ) , 
be = (2/1CCOS7, 2/icsin7, d 2 + c 2 — b 2 — a 2 — h 2 ) 



and 



G 



( 2ab + 2cd -2/iacos7 ^ 

2ab — 2cd —2hasm , y 
-2/1&COS7 -2/1&SU17 d 2 - a 2 - b 2 + c 2 + h 2 



(9) 



y— ZAiocos7 —zno sm 7 a — a — -\- c -\- a j 
Using these expression one can calculate five polynomial invariants and see that they are 
all even functions on 7. For instance, the three tangle is 

t = 4d^/(dh 2 - 4a6c) 2 + lQabcdh 2 cos 2 7. (10) 



Next one can inverse the obtained relations and express state parameters a, b, c, d, h and 7 
via polynomials \ip\ 2 , b A , &b> &c> £ and r(see equation (3.4) and comment below in Ref. [23]). 
But the inverse solution is not unique since both signs + or - are equally good for the 
gauge phase 7. Right here it appears a necessity to introduce the additional invariant i that 
distinguishes positive and negative values of the gauge phase 7. 

But to prove the theorem we do not need % at all. Indeed, from be = it follows that 
either h = or c = 0. Now i/j is a linear combination of four orthogonal product states and 
then the phase can be eliminated by appropriate local unitary transformations. Similarly, 
no additional invariant is needed if either = or b# = 0. 

A. Stationarity equations. 

The maximal product overlap g(ip) of a pure state can be expressed via the Bloch 
vectors and correlation matrix as follows [15] 

9= sup - [1 + x ■ b A + y ■ b B + x ■ (Gy)] , (11) 

x 2 =y 2 =l ^ 

where maximization runs over all unit vectors x and y. Note that x and y are the Bloch 
vectors of the local states \q A ) and \qs), respectively. 

By introducing Lagrange multipliers Ai and A2 that enforce unit vectors x and y one 
obtains the following stationarity equations: 

Gy + h A = \ lX , (12a) 

G T x + b B = X 2 y. (12b) 
Unknown Lagrange multipliers Ai and A2 are defined by the conditions 

|af = l, \y\ 2 = l. (13) 

These conditions are a pair of algebraic equations of degree six in Ai and A2. Evariste 
Galois's theory states that there is no general answer in terms of radicals and then the 
only possibility to find roots is the factorization of algebraic equation. The factorization 
of generic polynomials is impossible and therefore one cannot expect a simple closed form 
solution for arbitrary three-qubit states. Similar results have been obtained in Ref. [16], 
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where the authors derived a polynomial expression of degree 12 in the coefficients of a 
three-qubit pure state, the roots of which include the maximal product overlap of the state. 

Our main tasks are to solve the couple of equations (12) when the length of a some Bloch 
vector vanishes and show that g 2 = 1/2 in this case. It suffices to consider only the case 
be = since the analysis of the remaining cases is similar. Furthermore, from be = it 
follows that either h = or c = 0. These two cases will be analyzed separately in the next 
sections. 



B. Quadrilateral states 

In this section we consider the case h = 0. Now the correlation matrix G is diagonal and 
this fact essentially simplifies the stationarity equations. 

The state (7) is equivalent, up to local unitary transformations, to the state 

= a|100) + 6|010) + c|001) + (14) 

Fortunately this state is analyzed in detail in Ref. [15], where it is shown that an arbitrary 
entanglement measure is a fully symmetric function on state parameters a, b, c, d in this case. 
In particular, the three-tangle [21] of the state (14) is r = 16abcd and thus ip' is the GHZ 
state when a = b = c = d and is a W state when abed = [12] . 

The maximal product overlap of ip' is(up to factor 2) the circumradius of the cyclic 
quadrangle with the sides a, b, c, d and therefore this type of states can be categorized as 
quadrilateral states. The nearest product states of these states can be computed analytically 
by solving the stationarity equations (12). In particular, the local constituents of the nearest 
product state of if)' are [15] 

\ U ) = ^M±£SM, (is) 



\Qb) 
\qc) 



AS\/ad + be 
^Wd |0g) + Vrg^jlg) 

AS\/bd + ac 
y/rcTd\0c) + y/rJ'b |lg) 
ASy/cd + ab 
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where 



r a = a(b 2 + c 2 + rf 2 - a 2 ) + 26crf, 

r h = b(a 2 + c 2 + d 2 -b 2 ) + 2acrf, 

r c = c(b 2 + a 2 + d 2 - c 2 ) + 2abd, 

r d = d(b 2 + c 2 + a 2 - d 2 ) + 2abc 

and S is the area of the cyclic quadrangle a, b, c, d. 
From be = and h = it follows that 

c 2 + d 2 = a 2 + 6 2 . (16) 

Note that owing to this condition the correlation matrix G acquires a zero eigenvalue. But 
we were not forced to use it since the general solution for h = was already found in Ref.[15]. 
The condition (16) simplifies expressions (15) for the nearest product as follows 

Vbc\0 A ) + Vad\l A ) 
\Qa) = 7= ; = , (17) 

\Qb) 

\qc) 

The substitution of these expressions into Eq.(l) gives 

(c 2 + d 2 )ab + (a 2 + b 2 )cd 
9 — — / (lo) 
a/ (ad + be) (ac + bd) (ab + erf) 

Now from the identity 

(ac + bd)(bc + ad) = (c 2 + d 2 )ab + (a 2 + 6 2 )rfc 
and the normalization condition it follows that 

t = (19) 

C. The case h ^ 0. 

In this section we consider the case c = 0. Then the three-tangle is r = d 2 h 2 /A and the 
state is: the GHZ state when d = h and a = b = 0, a W state when h = and a biseparable 
state when rf = 0. 



v^IOa) + 


varf | 1a) 


yarf - 


h 6c 


v/aclOs) + 


v^IIb) 


Vac H 


-bd 


v^lOc) + 


Vab\l c ) 


Va6 H 


- cd 
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In general, if h 7^ then the system of stationarity equations are unsolvable since Eq.(13) 
yield generic algebraic equations of degree six. Moreover, the method developed in Ref.[16] 
also results a nonfactorizable characteristic polynomial of degree 12 for h ^ 0. In this reason 
the geometric measure of these states has not been investigated so far except particular 
cases [25-28]. 

Now we solve the stationarity equations explicitly when h ^ but be = 0. Fortunately 
when be = the correlation matrix G has a zero singular value and it is of crucial im- 
portance here. Owing to the existence of the zero singular value the sextic equations (13) 
can be factorized to the linear and quadratic equations. However, we use a singular value 
decomposition [29] instead to get rid of laborious algebra and find the answer as quick as 
possible. 



1. Singular value decomposition 



The requirements be = and c = impose the following condition 



d 2 = a 2 + b 2 + h 2 



(20) 



Furthermore, the gauge phase 7 can be eliminated by appropriate local unitary transforma- 
tions and we set 7 = for the simplicity. Then nonzero Bloch vectors are 



byi = ^^(sina, 0, cos a), b B = 6 B (sin/3, 0, cos/3), 



(21) 



where 



b A = 2aVh 2 + a 2 , b B = 2bVh 2 + b 2 , tana 



h 



The correlation matrix G is given by the simplified formula 



G 



( 



V 



2ab -2ha 
2ab 
-2Kb 2h 2 



The singular value decomposition of G is 



tan j3 



h 
b' 



(22) 



G = UDV + , 



(23) 
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where 



U 



cos a sin a 

1 
— sin a cos a 



and 





2/i 




, D = 


2ab 


, v = 




v y 





cos/3 sin/3 
1 



(24) 



H= ^{h? + a?){h? + b 2 ). 
It is easy to see that 

b A = b A Us, b B = b B Vs, (25) 
where s = (0, 0, 1) is the eigenvector of D with zero eigenvalue, i.e 

Ds = 0. (26) 
Now we define new unit vectors x' and y' as follows 

x' = Ux, y' = Vy. (27) 
Substituting expressions (23), (25) and (27) into stationarity equations (12) we obtain 

By' + b A s = Xix', (28a) 
Dx + b B s = X 2 y. (28b) 



2. Classification of the solutions 

Using equations (26) and (28) one can show that 

(D 2 - X^Dx' = 0, (D 2 - X^Dy' = 0. (29) 
Consequently there exist the following three types of solutions: 

1. The first type of the solutions are related to the zero eigenvalue of D and exist when 
both x' and y' are the zero mode of D. 

2. The second type of the solutions are related to the middle eigenvalue 2ab of D and 
exist when 

AiA 2 = (2ab) 2 . (30) 

3. The third type of the solutions are related to the largest eigenvalue 2/i of D and exist 
when 

AiA 2 = (2/x) 2 . (31) 
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3. Solutions of stationarity equations 



a. Zero mode solutions. From Eq.(29) it follows that stationarity equations (28) have 
the solutions Dx' = Dy' = 0. There are four these type of solutions: x' = ±s, Ai = 
±b A ,y' = is, A 2 = ±&b- In what follows we will omit all the solutions with negative 
Lagrange multipliers since only the solutions with Ai > and A2 > give a true local 
maximum. Then the zero mode solution of interest is 

x' = s, y = s, Ai = b A , A 2 = b B . (32) 

It gives the following product overlap 

gl = \{l + b A + b B ). (33) 

b. Nonphysical solutions. Consider now solutions given by Eq. (30). The scalar projec- 
tion of Eq.(28) onto s gives b A — X±(s ■ x') and b B = A 2 (s • y'), or 

b A b B = AiA 2 (s • x')(s ■ y'). (34) 

But b A b B > 4a 2 b 2 = AiA 2 , while (s • x')(s ■ y') < 1. Hence there are no physical solutions 
in this case. Sudbery et al. point out that some real roots of the characteristic polynomial 
have no associated singular vectors and does not mean a local maximum [16]. This is the 
case, the Bloch vectors of the local states are not real unit vectors and the maximal product 
overlap defined by (30) is not supported by a product state. It should be neglected. 

c. Main solution. Consider now solutions given by Eq.(31). We do not present the 
derivation of the solutions but describe the main steps. First, one parameterizes unknown 
vectors as follows: x' = (xx,X2,xz) and y' = (y±, 1/2, 2/3). Second, the condition AiA 2 = 4/z 2 
forces X2 = IJ2 = and factorizes the quartic equation into quadratic equations. 

We present only the solution with strictly positive Lagrange multipliers given by 



'6i + V 
b\ + 4/x 2 



A i = W it — a o = 2(a 2 + h 2 ), (35a) 



lb 2 B + 4// : 

and 



x > = 2 ^-sf^ = 2[il + h2) (35b) 
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Now one can put (36) into (28) and convince oneself that it is a solution of the stationarity 
equations with the Lagrange multipliers (35). 

The solution (36) gives the following product overlap 




D. The maximal product overlap. 

We computed the local maxima of the product overlap and now we would like to give a 
comment on specific points of entanglement parameters. The invariant be vanishes in the 
following two cases: either h = and c 2 + d 2 = a 2 + b 2 , or c = and d 2 = a 2 + b 2 + h 2 . It 
is easy to verify that the invariant t vanishes in both cases and and be are the left and 
right zero modes of G, respectively. On the other hand it is shown that all the five invariants 
are independent [19] and can be varied freely. What does this discrepancy mean? 

The explanation is that it can happen at the edge of the region of possible values of the 
invariants [30]. As a simple example, consider the three-dimensional set of points in the 
unit ball, x 2 + y 2 + z 2 < 1. The coordinates x, y, z are independent and can be varied freely 
in the interior of this region, yet we can prove that if z — 1 then x — y — 0. This is the 
case since be > 0. At minimal value of be the sextic invariant vanishes and the correlation 
matrix acquires a zero singular value. And it enables us to compute analytically maximums 
of the product overlap. 

The maximal product overlap for the case be = but h ^ should be defined as 

g 2 = max(# 2 ,# 2 ). (38) 

But 

b A + b B = 2aVa 2 + h 2 + 2bVb 2 + h 2 < 
a 2 + (Va 2 + h 2 ) 2 + b 2 + (Vb 2 + h 2 ) 2 = 1, 

therefore g\ < g\. Then 

g 2 = 1/2. (39) 

Hence we have shown that from be = it follows that g 2 = 1/2. Similarly, if either b A = 
or &b = 0, then g 2 = 1/2. The theorem is proved. 
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IV. FOUR- AND HIGHER-QUBIT STATES. 



An important problem is the generalization of the theorem to arbitrary n-qubit states. 
Unfortunately this cannot be done since the theorem ceases to be true when n > 4 [31]. 
Below we illustrate what is happening in this case. 

Consider first generalized GHZ states which can be written 

\GHZ) = cos 0|OO •••()) +sin0|ll---l) (40) 

in some product basis. The Bloch vectors are all equal in (40) and we denote them by 
b. Then g 2 (GHZ) = (1 + |b|)/2 [23] and thus the theorem is valid. By the way, using 
the Schmidt decomposition one can establish the same relation between g and b for pure 
two-qubit states. It means that pure two-qubit states with one single-qubit density matrix 
being completely mixed form a single orbit under local transformations and owing to this 
the theorem is correct. 

Consider now generalized W state which can be written 

\W n ) = ci|100...0) + c 2 |010...0) + ■ ■ ■ + c„|00...01) (41) 

in some product basis. It is not a trivial task to show that a zero Bloch vector forces g 2 = 1/2 
and viceversa in this case. The proof can be found in Ref. [20] and the theorem is true for 
generalized W states too. 

Consider now the four-qubit Dicke state given by 

\D) = 4=(|0011) + [0101) + [0011} + [1001) + [1010) + [1100}). (42) 
V6 

All Bloch vectors of this state are zero, but g 2 = 3/8^ 1/2 [14]. Hence the Dicke state is a 
counterexample showing that the straightforward generalization to higher-qubits is impos- 
sible. Perhaps the reason is the following. In the case of pure four-qubit states the full ring 
of local polynomial invariants is more complicated even in the presence of completely mixed 
states [32]. This means that in the best case the theorem can be somehow modified. 

V. DISCUSSION 

We have shown that all pure three-qubit states that possess a completely mixed one- 
particle density matrix have the same geometric measure of entanglement. This result raises 
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several questions. For instance, one may ask whether the theorem is valid for mixed states 
too. The answer is no. A simple counterexample is a completely mixed three-qubit state 
whose Bloch vectors are all zero but the state is separable. 

The most important question is whether the theorem can be extended to other entangle- 
ment measures suitable for arbitrary multipartite states. There are two sound arguments 
that the extension is indeed possible. First, the keystone idea of the proof is that at the 
edge of a quartic polynomial the sextic polynomial vanishes and the algebraically dependent 
invariant becomes redundant. And then two Bloch vectors become left and right zero modes 
of the correlation matrix. But this is a specific feature of local invariants of pure three-qubits 
and, therefore, is related directly to quantum entanglement. This feature has no relation 
to the geometric measure, or, more precisely, it is the same peculiarity for all entanglement 
measures. Hence any reliable entanglement measure should detect this peculiarity as a spe- 
cific point. Second, as it is pointed out in Ref. [15], all of the states with g 2 = 1/2, designated 
as shared quantum states, can be used as a quantum channel for the perfect teleportation 
and dense coding. Therefore they must possess the same amount of entanglement and thus 
a reliable entanglement measure should not vary on the manifold of shared quantum states. 
Again one can conclude that a completely mixed state should be a critical point for a good 
multipartite measure. 

Another important question is: what is happening with general n-qubit states? We have 
seen that both two- and three-qubit pure states have the same type of critical points that 
are the edge values of Bloch vectors. And when n > 4 the theorem ceases to be true. But 
our brief analyze in Sec. VI shows that the situation is more complicated. Indeed, for GHZ- 
and W-class states the theorem works well while it is wrong for Dicke states. This may 
indicate that there are different types of critical points at n — 4. And GHZ and W states 
have one type of critical points and Dicke states have another types of a such points. Then 
we would like to know how many types of critical points exist in the space of entanglement 
parameters of pure four-qubit states. And what is the role of those points or what kind of 
quantum phenomena are behind them. To clarify these points we need to analyze carefully 
the complete set of polynomial invariants of pure four-qubit [32] and two-qubit mixed [33] 
states. 

Finally we would like to discuss whether the inverse theorem is true. In the case of pure 
two-qubit states it is an easy task to show that if g 2 = 1/2 then b^ = bs = 0. But in the 
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case of pure three-qubit states the problem is open. The question is: does from g 2 = 1/2 it 
follow that bj^B^c — 0? Unfortunately we failed to prove or disprove the inverse theorem. 
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